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Abstract. We describe and discuss an experimental set-up which allows undergrad- 
uate and graduate students to view and study magnetic levitation on a type-I super- 
conductor. The demonstration can be repeated many times using one readily available 
25 liter liquid helium dewar. We study the equilibrium position of a magnet that levi- 
tates over a lead bowl inmiersed in a liquid hand- held helium cryostat. We combine the 
measurement of the position of the magnet with simple analytical calculations. This 
provides a vivid visualization of magnetic levitation from the balance between pure flux 
expulsion and gravitation. The experiment contrasts and illustrates the case of mag- 
netic levitation with high temperature type-II superconductors using liquid nitrogen, 
where levitation results from partial flux expulsion and vortex physics. 
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1. Introduction 

Meissner effect i.e., the expulsion of magnetic field from a superconducting material, 
has been extensively studied, both from the theoretical and experimental points of 
view [DElEllllElElElIHld} To explain analytically the observed results different 
methods are used [II1II2], which are, for some simple situations, tractable by students. 
In type-I superconductors, like lead or aluminum, the Meissner effect provokes a total 
expulsion of the magnetic flux, and a magnet located on top of one of these materials is 
rejected due to repulsion. Nevertheless, typical demonstrations of magnetic levitation 
on a superconductor show a magnet levitating in apparent equilibrium on top of a 
superconductor. Of course, this occurs when using cuprate high Tq superconductors, 
which are type II materials, where the magnetic field enters in the form of quantized 
vortices. When a magnet is brought close to a type-II superconductor, the behavior 
observed is the result of Meissner flux expulsion combined with vortex induced attraction 
or expulsion. Moreover, the magnetic history of the superconductor, with possible flux 
trapped, significantly influences the position of the magnet. Many related videos can 
be found in the web (see, for instance, references [131 [H [ISl HSl [T7]). 

Vortices are pinned to imperfections in the sample (vacancies, impurities, 
dislocations, etc) and they tend to be arranged in a regular lattice (see, for instance, 
[181 [TSj and references therein). When the vertical position of the levitating magnet 
is altered just by pushing or unfastening, the vortex distribution is changed in the 
superconductor. This can occur as long as the pinning is not very strong. The same 
applies to lateral motion. In this way, in Figure [l|^a) we can see a magnet located just 
above the center of the superconductor, with a symmetric distribution of vortices. In 
Figure [T]^b) we push the magnet towards one side of the superconductor, and we obtain 
again an equilibrium position, although the vortex distribution has changed (we neglect 
demagnetizing effects due to superconductor's edges). On the contrary, when pinning 
is really strong the magnet can be difficult to move, and any small displacement is 
corrected as it tends to return to its initial position. Lateral stability can be taken for 
granted when using type-II superconductors, which is not the case for type-I, due to 
repulsion and the lack of flux penetration below the critical field, unless the intermediate 
state appears. This regime can arise in type-I superconductors, when the geometry of 
the superconducting body provokes that the density of flux lines is not homogeneous 
around its surface (for instance, for a sphere the density of flux lines would be higher 
at the equator and zero at the poles). Hence, the local magnetic field can be above 
the critical value, Hq, in some parts, whereas it remains below in other regions around 
the sample. This leads to the formation of normal and superconducting domains, whose 
frontiers are always parallel to the applied field, H. On the contrary, in the cross-section 

I There are plenty of books where superconductivity phenomena are explained in depth, for instance 
in reference [3]. A visual demonstration of superconductivity phenomena performed at the Michigan 
State University can be found in a set of 6 videos. The address for the first one is given in reference 
[TU] , Meissner effect is described in segments 3 and 4. 
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perpendicular to H, the distribution of normal areas exhibits peculiar and irregular 
domain patterns, depending on the geometry (sphere, foil, thin wire, etc). Note that 
the stable coexistence of these domains requires that the field in the normal regions is 
just Hq. Below this value the superconductivity would be regained, and above it the 
normal zones would propagate towards the superconducting adjacent regions [§j 

In Figure [l](c) we can see a magnet levitating above the center of a type-I 
superconductor. Vortices do not exist because all the magnetic flux is expelled, and so 
there is pure repulsion. The lateral stability depends on the particular geometry of the 
system and the magnet can even exhibit a precession around a rotation axis [21 [3] . This 
is avoided if the type-I superconductor is fabricated with a concave shape, like a bowl. 
The horizontal components of the radial repulsive forces are directed to the magnet in 
such a way that it can be stabilized around the center of the bowl [2S]. In Figure [l|^d) 
we have depicted a scheme of a magnet levitating over a type-I superconducting vessel. 
The horizontal components of the repulsive force, F^, keep the magnet in an equilibrium 
position at the center of the vessel. Note that the absence of a vortex lattice implies 
that there is only one equilibrium position. 

In this work, we present a simple set-up to view magnetic levitation with a type- 
I superconductor. We propose an analytical procedure to calculate the equilibrium 
position of a cylindrical magnet that levitates over a vessel made from lead, whose 
critical temperature is 7.2 K. We provide simple calculations, using basic magnetostatics, 
to explain the observed stable position of the magnet on top of the superconductor. 



2. Experimental procedure 

As mentioned above, the first step was the realization of the experiment with a vessel 
made from lead and a cylinder made from a NdFeB alloy (see the complete videos, 
recorded during the realization of the experiment, at [261 122] )• The vessel was inside a 
glass hand-held cryostat that has a double wall with vacuum insulation ^ Figure [2|a) 
shows a photograph of this dewar and the lead bowl. 

As the critical temperature of lead is 7.2 K, liquid helium is needed. Normally, this 
cryogenic liquid must be handled with care, and protective gloves and goggles must be 
worn. Figure |2](b) shows how to hold the cryostat when pouring liquid helium. The 
siphon is of a common type and can be supplied with the helium bottle (safe handling 
and storage of liquid helium is briefly described in [2^). In addition, the latent heat of 
evaporation of helium-4 is much lower than that of liquid nitrogen (20.6 J/g and 199 
J/g, respectively, at their boiling temperatures and so it evaporates very quickly in 
a simple dewar as the one used here. 

§ Further insight into this subject and nice images of the intermediate state can be found in 

[211 Ea EH El- 

II More information about how to obtain small glass dewars where liquid helium is preserved for a few 
minutes can be obtained from authors. Made to measure glass tubes can be purchased at [128 . 
% A comprehensive summary of the aspects involved in the handling of liquid helium can be found in 
the 5th chapter of reference [30] , 
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Magnet under 
repulsive forces 




Figure 1. (a) Magnet levitating over a type-II superconductor. Pinned flux lines 
keep the magnet in a stable position, (b) Once the magnet is forced to move, another 
equilibrium position can be found, (c) Magnet over a type-I superconductor. As the 
magnetic flux is completely expelled from the superconductor, the magnet is not in 
a real stable position, and it may rotate or oscillate, (d) Magnet over a vessel made 
from a type-I superconductor. The horizontal components of the radial force allow a 
certain lateral stability, as depicted in the top view on the right. 



We tried two different variations of this experiment. In botfi cases, the lead vessel 
was located in the bottom of the hand-held cryostat and a glass tube was used to guide 
the magnet towards the vessel. Nevertheless, some aspects were different, these being: 
(i) the magnet was let to slide down through the tube and towards the centre of the 
vessel after the cooling of the whole set-up. This is the experiment to be compared 
with our analytical model; (ii) a less powerful magnet was released through the tube 
but away from the symmetry axis of the vessel. Let us see these cases separately: 

(i) We first poured some liquid helium inside the cryostat. After some seconds, boil 
off relaxed and the lead bowl reached the full superconducting state. Hence, the 
glass tube was introduced to guide the magnet to the symmetry axis of the vessel 
and keep it in a vertical position. Once the magnet was let to slide inside the tube, 
it attained the lead bowl but, immediately, bounced and remained stable at an 
equilibrium position, about 1 cm above the brim of the vessel (as can be seen in 
Figure [Sj^a), using as reference the height of the lead bowl, equal to 2.6 cm). As 
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Figure 2. (a) Photograph of the hand-held cryostat we used for the reahzation of the 
experiment. (b) The cryostat can be handled easily when pouring liquid helium, but 
using protective gloves is strongly recommended. 

could be expected, displacing the glass tube did not alter the equilibrium distance. 
As noted above, the lack of pinning provoked that there was just one stable position, 
which could not be altered. In the photographs it can observed that the magnet 
was slightly tilted, probably due to an inhomogeneous field distribution. This is not 
surprising, taking into account that the vessel was not totally regular, and hence 
the horizontal components of the repulsive forces could be somewhat unbalanced. 
After about 30 seconds, once the hquid helium was completely evaporated, the lead 
transited to the normal state, and the magnet fell down the tube and inside the 
vessel, as shown in Figure |3](b). 

(ii) In this experiment a small magnet was thrown in the cryostat, away from the 
symmetry axis of the vessel. In this case, it was completely repelled, and it was 
leaning against the cryostat walls even when this was rotated. This was the expected 
behaviour, taking into account the lack of symmetry of the forces when the magnet 
is apart from the vessel's centre. 

Concerning the refrigeration costs of these experiments, we can make a rough 




Figure 3. (a) Magnet levitating in an equilibrium position over the lead vessel, (b) 
Magnet falling down once the bowl leaves the superconducting state. 
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estimation of the liquid helium wasted during the whole run. First of all, the volume 
of liquid required to cool the magnet (the biggest one) from 300 K to 4.2 K can be 
calculated just by using its latent heat. Assuming that the Debye temperature of the 
magnet is Gj^ ~ 300 K, and taking into account the Debye tables (see [30], section 
11.2.5 and Table B2), the energy to be removed from the magnet is roughly 313 J. As 
the latent heat of liquid helium is 2.6 kJ/1, about 120 ml of helium would be needed 
to cool the magnet. Nevertheless, in a real great deal of the cooling process is 

undertaken by the enthalpy of the gas (about 200 kJ/1 between 4.2 K and 300 K) and 
a lot of liquid is saved [31J. In fact, we estimate that we lose less than 30 ml of liquid 
helium each time we cool down the vessel, the glass duct and the magnet. 

3. Calculation of the equilibrium position by using basic electrodynamics 

3.1. Analytical model 

We present in this section a simple calculation to estimate the equilibrium position of 
the magnet on top of the superconductor. This calculation can be used in e.g., an 
electrodynamics course, together with the demonstration of the experiment. 

As a first step, let us consider a cross-section of the system under study and a 
proper set of distances and dimensions, as depicted in Figure |4j In this scheme, R and 
L are the magnet's radius and length, respectively, h the height of the lead vessel, a 
its inner radius, b the outer one, t the thickness of its base and, finally, zi the distance 
between the bottom of the vessel and the magnet. The origin of coordinates is taken 
where the z axis intersects the bottom of the vessel. The numerical values of all the 
parameters involved in the calculation of zi are gathered in Table [1} 

Table 1. Dimensions of the magnet and the lead vessel. 

Parameter and units Numerical value 

Magnet radius, R (mm) 3 

Magnet length, L (mm) 15 

Magnet mass, m (g) 3 

Vessel inner radius, a (mm) 13 

Vessel outer radius, b (mm) 15 

Vessel height, h (mm) 26 

Vessel base thickness, t (mm) 2 



In order to proceed with the calculations, we are going to consider a magnet having 
a constant magnetization oriented along the z positive axis, i.e. M = Mu^ Q The 

+ To avoid carrying unnecessary subscripts, we are not going to use any special distinction for the 
magnitudes related to the magnet. The ambiguity is broken by the subscript "vessel" when referring 
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= 0, z = 0) Lead vessel 



Figure 4. Cross-section of tlie system magnet-vessel and definition of the dimensions 
and distances which are going to be relevant in our calculations. The origin of the 
coordinate system is also indicated. 



magnetization in the vessel can be calculated from the magnetic induction of the magnet 
as 

B 

and it is related to the surface magnetizing current density, J^^^g, through 



M 



vessel 



(1) 



JmS = n ^ ^vessel- 



(2) 



Note that in our superconducting vessel the thickness of this thin layer of current 
is the London penetration depth, and also that the constant magnetization precludes 
the existence of any volumetric current density i.e., Jj^y = V x M^ggg^^ = 0. 

The repulsion force between the magnet (with constant induction B) and the vessel 
can be expressed as 



Fm = I^JjnS 



X BdA, 



which, in turn, must balance the force exerted by gravity, so 

Fm = """^magnetS- 
to the lead bowl. 



(3) 



(4) 
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In order to solve Equation ([3]), we need first to find a convenient expression for 
the magnetic induction. Taking into account that there are no free currents inside 
the magnet, the magnetic induction can be expressed as the gradient of a magnetic 
potential i.e., B = — /UoV0m- Outside the magnet, this potential can be calculated by 
using Equation ([S]) [32j: 

<Pm{r) = ^ f?^dA', (5) 



An Js ^ 

where it has been used the magnetic pole surface density, defined as crm(i"') = n ■ M(r'). 
Due to the magnetization orientation, am is just defined on both the top and bottom 
faces of the cylinder. Therefore, we will have = Uz ■ M = M and o"]^^^'^^ = 

— Uz ■ M = —M. We have not included the term corresponding to the magnetic 
pole volume density, pm(i"') = —V ■ M(r'), which is zero because of the constant 
magnetization. It is to be remarked that this approach, based on magnetic pole densities, 
is useful just outside the magnet, and that we would get erroneous results by applying 
Equation ^ within its own volume. 

Then, the scalar magnetic potential can be written as: 

M fR p'dp' 



M fR p'dp' 
2 ^0 ^{p-p'y + iz,-zf 

where we have already performed the integration over the angular coordinate {J^^ dip' = 
27t). Note that the integral is calculated over the source radial coordinate, p', which 
varies along the cylinder radius, whereas z' is fixed and equal to Zi for the bottom 
cylinder's base and Zi + L for the top one. 



By using a Taylor series and some convenient substitutions (see Appendix A.l), the 
scalar magnetic potential can be written in the following way: 




{Zi - 2)3 (zi+L- Z] 



3 



(t - 3^ ^+ j - T (., + L-.)(.,-.) - 

Now we can find out the magnetic field density by deriving this magnetic potential. 
As we are using a cylindrical coordinates system, the gradient components will be: 

\ dp ^ p dip ^ dz J ' 

The magnetic potential has not any angular dependence in this problem, so there are 
only two components of the magnetic field density, these being 

Bp = ~/^o-^, and = -/^o-^- (9) 
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1 



(zi + L - z 



|3 



(10) 



The other component, B^, will be not calculated now because we will not use it to obtain 
the force in the Uz direction, as we will see below. 

Now that we have a suitable expression for B, we can calculate an analytical 
expression for the repulsion force. Combining Equation ^ and Equation (|3]), this 



force can be written as an integral over the whole vessel surface (see Appendix A. 2 for 
details) : 



m 

/io Js 



in Js 



fB-n)B)dA. fill 



This equation can be expressed as a sum of integrals corresponding to every face of the 
vessel. Then, taking into account the values of the magnetic field induction components. 



the repulsion force in the Uz direction is (Appendix A.2): 

2n 



Fm = -Uz— / pBl{p,z = 0)dp- 
Po Jo ^ 

27r 

Uz — / hBp{p = b, z)B^{p = b, z) dz + 
Po Jo 

Uz — / pB^p{p, z = h) dp + 

Pq J a 

Uz — / CiBp{p = a, z)Bz{p = a, z) dz + 

Pq Jt 

Uz- / pB%p,z = t)dp. (12) 
Pq Jo ^ 

The first member corresponds to the outer base, the second to the outer lateral face, 
the third to the brim, the fourth to the inner lateral face and the fifth to the inner base. 
Notice that we have made explicit the value for every non-integration variable. 



The analytical solution of Equation (12) does not yield a compact and easy-to- 



handle result, so we propose an approximation, this being that the magnetic field 
scarcely varies within the small thickness of the walls. In this way, it is assumed that 



the second and fourth terms in Equation (12) are virtually equal, except for their sign 



and so its sum is close to zero, and the same applies to the first and the last terms. This 
is the reason why we do not need any analytical expression for the component -B^, as 
stated above (later on, we will analyse the quality of this approximation). Finally, only 

* It must be remarked that this value of Bp is vahd just in the vicinity of the magnet, as commented 
in Appendix A.l 
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the third term of this equation survives, and we have a repulsion force given by (see 



Appendix A.3): 



Fm~/iOoM2j-(i?,a,6) 

^ , (13) 



1 V 



J-'(i?, a, h) being a simple function of R, a and h. 
3.2. Numerical result 

The magnetic induction created by the magnet is 5 = 4400 Gauss [f} and so M 



-B(/io)^^ = 0.35 MA/m. Finally, by making Equation (13) equal to the force exerted 
by gravity (Equation (|4])), and by using the numerical values of Table [l| we get the 
following relation: 



1 



^ =1.5-10^^ (14) 

V(^i- 0.026)3 (^1- 0.011) V ^ ' 

Solving this equation by simple iteration yields a result of Zi = 3.52 cm. Taking into 
account the definition of zi and h given in Figure |4], we find that the distance between 
the brim of the vessel and the bottom of the magnet is — /i = 0.92 cm. Despite 
the magnet did not levitate exactly at the center of the vessel and it was somewhat 
tilted, this result is in excellent agreement with the experimental value, as can be seen 
in Figure |3](a). 

3. 3. Repulsion force calculated without approximations 



For convenience, the exact solution of Equation (12) has been carried out by numerical 
integration. Figure [5] shows the comparison of the force obtained just with the 
third integral (dashed line) and with the whole equation (solid line), as well as their 



intersections with the force exerted by gravity (dotted line) (see Appendix A. 4 for the 
analytical expression of B^). It can be observed that the curves do not coincide, but 
there is a 2 mm difference between both crossing points. Hence the new equilibrium 
distance would he zi — h = 0.71 cm. 

The reason for the lack of accuracy of our approximation is that the integrals 
containing the component cannot be neglected (the first and fifth integrals do yield 
contributions which are several orders of magnitude smaller). Despite this error, the 
results are not so much different, and our calculation, intended to be as easy-to-handle 
as possible, allows obtaining quite a reasonable approach to the equilibrium distance. 

tt The characteristics and prices of these items can be found in [33j. 
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Figure 5. Repulsion force calculated just with the third integral of Equation (12) 
(dashed line) and with the whole set of terms (solid line). They yield intersections 
with the gravity force (dotted line) which differ in about 2 mm. 



3.4- Considerations on the radial force 

Provided that the magnet is located along the symmetry axis of the vessel, the resultant 



of the radial force is zero (see appendix Appendix A. 2 to see how this is derived 



from calculations). This was already commented in the introduction and Figure [T|d), 
regarding the lateral stability of the magnet during the levitation. Nevertheless, if the 
magnet is out from the symmetry axis the scenario is different. In that case, there would 
be a non-zero resultant radial force, because of the unbalancing of the individual radial 
components, and the magnet would be repelled out from the vessel. This effect would 
explain the second case considered in section |2| where a small magnet was expelled from 
the center of the vessel towards their borders. 

Despite this, we observed that the cylindrical magnet we used for the experiment 
involved in the calculations was not repelled when its position was not exactly at the 
center of the vessel. This can be explained considering that it is a handmade bowl and 
not regular in shape, so the real equilibrium point for the magnet can be expected to 
be shifted to one side, or there may be different places where the magnet is more or less 
stable. In addition, the fact that the magnet is somewhat tilted with respect to the z 
axis provokes changes in the balance of forces |34j . 

On the other hand, it is possible that the off-axis equilibrium point is somehow 
related to the intermediate state, already mentioned in the introduction. This could 
be reasonable, taking into account that the cylindrical magnet is much stronger than 
the small magnet we used when we observed the total repulsion against the glass tube. 
In that scenario, some magnetic flux could have penetrated into the superconducting 
vessel, yielding a somewhat unpredictable equilibrium position. 
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4. Conclusions 



We have reproduced and recorded in our laboratory the levitation of a magnet 
over a type-1 superconductor by using a simple hand-held cryostat, and we have 
obtained from simple electromagnetism the analytical equation for the equilibrium 
position. Experiment provides for a visual discussion of the always fascinating levitation 
phenomenon. The pedagogical ingredients complement the popular levitation with 
type-II high Tc superconductors. Our experiment includes thermodynamics (properties 
and handling of cryogenic liquids) and tractable electrodynamics. Handling of liquid 
helium is greatly simplified with a hand- held small sized cryostat. The experiment 
can be carried out without a great cost in liquid helium. Using 25 liters of liquid 
helium (approximately 250 euros), we have shown the Meissner effect and discussed the 
properties of liquid helium to more than 200 secondary school students each morning 
during 3 days. This experiment also serves as an introduction to advanced solid state 
concepts (superconductivity), if combined with magnetic levitation in high Tq materials. 
It can be used to discuss the applications of phenomena which are enabling industrial 
solutions of great importance to efficient and sustainable energy handling, such as 
cryogenics and superconductivity. 



Appendix A. Additional details on the calculations 

In this appendix we present some details on the different steps that are not strictly 
necessary to follow the calculation. 



Appendix A.l. 

Equation ([T]) was obtained by applying a Taylor series around x = 

' ' ^ (A.i: 



\Jx + a ^Ja 2a^/2 

and taking into account that x = {p — p')^, and a = {zi + L — zY for the first term of 
Equation ^ and a = {zi — zY for the second on^ff] In this way, we have Equation 



'^-"Ti(^r3^-(., + L-.)3jio P^P-P^ 

2 {zi + L — z)[zi — z) Jo 
from which Equation ([T]) results after integration. 

ft Note that this approximation impHes that p ^ p' i.e, the field and source points are very close. 
Despite its apparent roughness, it yields good results for distances about the dimensions of our set-up. 
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Appendix A. 2. 



Equation (11) was derived by substituting 3ms = — — n x B into Equation (3), and 



Mo 



then applying the vectorial identity (a x b) x c = (c ■ a)b — (c ■ b)a. Then, we expanded 
this last equation in terms of the different contributions due to each face of the vessel: 



1 /'27r rb 

/ / (^^(-u,) - (B ■ {-%)B))pdpdip + 
Jo Jo 



m 



1^0 



1 .2^ [^(B'^vL - (B ■ Vi,,)B)hdzdip + 
Jo Jo 

1 f " [\b^u^ - (B ■ u^)B) pdpdifi + 

Jo Ja 



1^0 

I r2TT i-h 



r j (52(-u ) - (B- (-u„)B))adzd(^ + 
Pq Jo Jt 

— r Hb^u, - (B ■ u,)B)pdpdip. 
Po Jo Jo 



(A.3) 



In addition, it was taken into account that = Bf, + B'f, B ■ 



B • Uz = Bz, as well as Jq^ u^dip = 27rUz, and J^^ Updip = (let us remember that 
Up = cos(/)Ux + sin<y9Uy, which yields zero when integrated between and 27i). 

Appendix A.3. 



Equation (13) was obtained when substituting the calculated expression for Bp into the 



third term of Equation (12), as commented, yielding 



TT 



Fm ^ f^o^M' 



Xzi-hf (zi + L-hf 



pi^-R'-R'pj dp. 



After integration, it was obtained Equation (A. 5) 



(A.4) 



TT 



m 



9 9 4 



(A.5) 



Notice that the term involving powers of -R, a and 6, is the function a, h) of Equation 
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Appendix A.4- 

The component B^, resulting from direct derivation of Equation ([T]) with respect to z, 



IS 



B. 



3M 



4 



— P H P" 

3 ^ 2 ^ 



,zi - z 



1 

[zi + L 



M 



+fio—LR' 



{zi - z){zi + L - zf {zi- zf{zi + L- z) 



(A.6) 



Although calculating the second and fourth integrals of Equation (12) is not 



difficult, the process is tedious and the resulting analytical equation is very long and 
not easy to handle. For that, as a first approximation, only the third integral was used 
to obtain the repulsion force. 
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